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Abstract. We study the backreaction on the mean field geometry due to a non- 
conformal quantum field in a Robertson- Walker background. In the regime of small 
mass and small deviation from conformal coupling, we compute perturbatively the 
expectation value of the stress tensor of the field for a variety of vacuum states, and 
use it to obtain explicitly the semiclassical gravity solutions for isotropic perturbations 
around de Sitter spacetime, which is found to be stable. Our results show clearly 
the crucial role of the non-local terms that appear in the effective action: they cancel 
the contribution from local terms proportional to the logarithm of the scale factor 
^ , which would otherwise become dominant at late times and prevent the existence of 

a stable self-consistent de Sitter solution. Finally, the opposite regime of a strongly 
\Q . non-conformal field with a large mass is also considered. 

(N 

o 

^ ' 1. Introduction 



Our present understanding of cosmology assumes that the universe underwent a short 
period of accelerated expansion known as inflation [1, 2, 3, 4, 5]. The inflationary 
scenario has been remarkably successful in explaining the observed anisotropies of the 
cosmic microwave background [6, 7, 8, 9]. In most inflationary models the accelerated 
expansion phase is close to but never exactly de Sitter and this phase eventually ends 
when the kinetic energy of the inflaton field driving inflation starts to dominate over 
the potential term. On the other hand, observations of distant supernovae indicate that 
the universe is presently undergoing a period of accelerated expansion [10, 11] that may 
be driven by a small non- vanishing cosmological constant [12, 13, 14, 15]. If that is 
the case, the geometry of our universe would tend to that of de Sitter spacetime at 
sufficiently late times. Thus, a detailed knowledge of the physics associated with de 
Sitter space may play a key role in understanding both the very early universe as well 
as its ultimate fate. Furthermore, it is conceivable that studying a possible screening of 
the cosmological constant driving de Sitter space, due to quantum effects, could shed 
some light on the huge fine-tuning problem that the current value of the cosmological 
constant seems to pose. 
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An open question which has recently received increasing attention is whether the 
quantum fluctuations of the metric and the matter fields in de Sitter space can give 
rise to large backreaction effects on the mean background geometry. It has been 
argued that in pure gravity with a cosmological constant the infrared effects due to two 
graviton loops and higher-order radiative corrections could lead to a secular screening 
of the cosmological constant [16, 17]. There have also been proposals that a significant 
screening of the cosmological constant could appear in chaotic inflationary models at one 
loop when both the metric and inflaton field fluctuations are considered [18, 19, 20, 21]. 
In all these cases the quantum fluctuations of the metric play an essential role. However, 
whenever the metric perturbations are quantized, one needs to confront the problem 
of defining proper diffeomorphism-invariant observables in quantum gravity [22], even 
when treated as a low-energy effective field theory. In particular one needs to make 
sure that the secular screening found in the analysis mentioned above is not simply 
a gauge artifact. As a matter of fact, it was shown in Refs. [23, 24] that when a 
suitable gauge-invariant measure of the expansion rate was considered, the screening 
effect previously found in chaotic inflationary models was not actually present (at least 
for single field models). Similarly, a recent reanalysis of the pure gravity case which made 
use of a diffeomorphism-invariant measure of the change of the expansion rate indicated 
the absence of secular effects to all orders in perturbation theory [25] (although this 
conclusion is still subject to certain debate [26]). 

In fact, the study of backreaction effects from quantum matter fields in de Sitter 
spacetime has a long history. Fischetti, Hartle and Hu employed effective action methods 
to study the quantum backreaction of conformal fields on the dynamics of a Robertson- 
Walker (RW) spacetime [27]. Shortly afterwards Starobinsky showed that the vacuum 
polarization effects of a large number of conformal fields could drive a de Sitter expansion 
stage without the need for a classical cosmological constant and that this de Sitter 
solution was unstable under RW-type perturbations (for a certain sign of one of the 
parameters in the effective action), which provided a "graceful exit" mechanism towards 
a standard cosmological evolution [28, 29]. This scenario is often known as Starobinsky 
(or trace anomaly) inflation. However, it has later been argued that within an effective 
field theory (EFT) approach (which provides a natural framework for understanding 
this kind of calculations) the higher-order curvature terms in the effective action should 
be treated perturbatively [30, 31, 32, 33]. When doing so, the de Sitter solution entirely 
driven by the vacuum polarization of the conformal fields is no longer present [34], f 
and a de Sitter solution driven by a cosmological constant A is stable even when the 
vacuum polarization effects from the conformal fields are included. More recently, it 

X One could argue that since the Hubble radius in Starobinsky inflation is of order >/Nlp (where N is 
the number of conformal fields), for a sufficiently large N the quantum gravity corrections should 
be small. Nevertheless, the same kind of argument would imply the existence of instabilities for 
perturbations around flat space with a characteristic time-scale of order vf/p [32]. Moreover, in 
particular implementations of Starobinsky inflation where the AdS/CFT correspondence can be applied 
[35] , the scale where higher order corrections become relevant and the low-energy effective field theory 
expansion breaks down is actually y/Nlp rather than lp. 
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has been argued that the backreaction from non-conformal fields can lead to significant 
deviations from a de Sitter solution (in some cases this has been analyzed in the context 
of trace anomaly inflation [36, 37], but it would also apply to a A-driven de Sitter 
spacetime) . Similarly, it has also been concluded that graviton one- loop effects can lead 
to significant deviations from de Sitter for the background dynamics [38, 39]. Since 
the metric perturbations are quantized in this case, there are gauge ambiguities, as 
mentioned above, which need to be understood and properly addressed. However, if 
one temporarily ignores this point, the resulting effective action which governs the 
background dynamics has the same form as that describing the vacuum polarization 
effects of non-conformal quantum matter fields. 

Here we reanalyze this problem by explicitly solving the backreaction on the mean 
gravitational field due to the quantum effects of a non-conformal scalar field when 
the quantum fluctuations of the metric are not considered. This kind of one-loop 
calculation is entirely equivalent to studying the corresponding backreaction problem 
in the semiclassical gravity framework [40, 41, 32] by solving the semiclassical Einstein 
equation, which includes the suitably renormalized quantum expectation value of the 
stress tensor operator acting as a source. Specifically, in our calculation we assume the 
presence of a cosmo logical constant, which would lead to a de Sitter solution in the 
absence of quantum effects, and simplify the problem by focusing on RW geometries, 
corresponding to spatially homogenous and isotropic states of the quantum field. 

There exist relevant antecedents to our analysis in the context of both quantum 
field theory in a fixed curved spacetime (when the backreaction of the quantum fields 
on the spacetime geometry is not taken into account) and semiclassical gravity. The so- 
called Bunch-Davies vacuum [42, 40] for fields in de Sitter is a state invariant under all 
the isometries of de Sitter space, which is maximally symmetric. The renormalized 
expectation value of the stress tensor operator for that state is proportional (with 
a constant factor) to the metric and, therefore, its contribution to the semiclassical 
Einstein equation has the same form as a cosmological constant term, which allows 
the existence of self-consistent semiclassical de Sitter solutions [43, 28, 44]. More 
importantly, it was shown in Ref. [45] that for fields with a wide range of mass and 
curvature-coupling parameters evolving in a given de Sitter spacetime, the expectation 
value of the stress tensor for any reasonable initial state tended at late times to the 
expectation value for the Bunch-Davies vacuum, where by reasonable states one means 
states with the same ultraviolet behavior as the Minkowski vacuum, i.e. with essentially 
no excitations at arbitrarily high frequencies (technically they are known as fourth- 
order adiabatic states [40, 46]). This result can be intuitively understood as follows: 
the exponential expansion will redshift any finite frequency excitations of the Bunch- 
Davies vacuum so that their contribution to the stress tensor will tend to zero at late 
times. This result suggests that even when taking into account backreaction effects, 
perturbations around de Sitter will be redshifted away and at late times the spacetime 
geometry will approach de Sitter space with an effective cosmological constant which 
includes the contribution from the expectation value of the stress tensor for the Bunch- 
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Davies vacuum. Such an expectation is supported by the analysis of [47, 48, 49], where 
the linearized semiclassical Einstein equation for spatially isotropic perturbations was 
considered. The asymptotic behavior of its solutions at late times was analyzed and 
found to be stable. 

Our approach enables us to obtain explicit and relatively simple results for the 
solutions of the semiclassical backreaction equation at all times, rather than just study 
its asymptotic stability. This can be achieved for the two opposite regimes of weakly and 
strongly non-conformal fields. The weakly non-conformal case corresponds to fields with 
a mass much smaller than the Hubble parameter and a small deviation from conformal 
coupling to the curvature. In that case one can treat perturbatively the mass and the 
parameter characterizing the deviation from conformal coupling. Our results for the 
semiclassical solutions, given by (67) and (79), confirm the stability of de Sitter for 
weakly non-conformal fields. This extends our results in an earlier study restricted to 
massless fields with a slightly non-conformal coupling to the curvature [50], where a 
similar method was employed. On the other hand, for the strongly non-conformal case, 
corresponding to a Compton wavelength much smaller than the typical curvature radius, 
one can introduce a quasi-local approximation that gives rise to a local expansion for the 
effective action involving positive powers of the curvature divided by the mass squared. 
We find that through order 1/m 2 no non-trivial spatially isotropic perturbations are 
allowed. These approximation schemes do not provide information on the intermediate 
regime of masses comparable to the Hubble parameter, but make it possible to obtain 
explicit expressions with a simple form for the other two regimes. One of the advantages 
of our explicit results is that they clearly show the crucial role played by the non- 
local terms (particularly for the weakly non-conformal case). Indeed, the effective 
action and the backreaction equation that one can derive from it exhibit local terms 
proportional to the logarithm of the scale factor which become dominant at late times 
and can cause substantial effects. However, these terms are canceled out by a similar 
contribution from the non-local terms. This is an important point which seems to have 
been unnoticed in previous studies and can have significant implications on the validity 
of local approximations where the non-local terms are neglected. 

The plan of the paper is the following. In section 2 we introduce the model 
and the class of initial states that we will be considering. We also explain how 
to obtain the expectation value of the stress tensor from the influence action and 
how to renormalize the divergences that appear in the influence action. Finally, we 
derive the backreaction equation and its solutions in terms of the expectation value 
of the stress tensor. The effective action for weakly non-conformal fields is computed 
perturbatively in section 3 and explicit results for the stress tensor and the solutions 
of the backreaction equation are obtained for massless non-conformally coupled fields 
and massive conformally coupled ones. A more detailed exposition of certain technical 
aspects can be found in [50], where the massless non-conformally coupled case was 
studied. In section 4 we describe how to obtain the influence action for strongly non- 
conformal fields with a large mass using a quasi-local expansion, and their backreaction 
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effects on small perturbations around de Sitter. We also explain the relevance of these 
results in order to compare with certain dark energy models based on the renormalization 
group (RG) running of the cosmological constant. Finally, we discuss the implications 
of our results in section 5. The relationship between the class of states that we consider 
and the fourth-order adiabatic vacua for RW spacetimes is explained in the appendix. 
Throughout the paper we use natural units with h—c—1 and the (+, +, +) convention 
of [51]. 

2. General model 

We consider a spatially flat RW metric, which in conformal time has the form 

g^u = a 2 {v)v^, (1) 

where r]^ = diag(— 1, 1, . . . , 1) is the n-dimensional Minkowski metric (we use arbitrary 
dimensions for the moment in order to perform dimensional regularization later on), and 
a(i]) is the scale factor in terms of the conformal time 77. We will study the dynamics 
of a free quantum scalar field of mass m in such a spacetime, whose action is 

S[a, ®] = -\J d Hxan KV^^ + [m 2 + + v)R}<$> 2 } , (2) 

where the dimensionless parameters £ c = (n — 2)/[4(n — 1)] (equal to 1/6 in four 
dimensions) and v give the coupling to R, which is the Ricci curvature scalar associated 
with the metric (1). It is given by 

„ / n ( H n — 4 d 2 \ 

fl = 2 (»- 1 »U+— ?)■ (3) 

Here and throughout the rest of the paper overdots denote derivatives with respect to the 
conformal time, i.e. ' = d/drj. The case v = is the so-called conformal coupling. If the 
field is massless and conformally coupled, then it is said to be conformally invariant, or 
simply conformal. In this paper we will be interested in the general case of non-conformal 
fields. It is convenient to introduce the rescaled scalar field <fi(x) = a^ n ~ 2 ^ 2 (7])^(x). In 
terms of this rescaled field, the action simplifies to 

S[a, <t>] = -\J dUx {v^d^dA + M 2 2 ) , (4) 

where M 2 is a function of the conformal time given by M 2 (r]) = a 2 (r])[m 2 + vR(rj)\. 
Setting the variation of this action with respect to equal to zero yields the dynamical 
equation for the rescaled field, 

(rrd^ - M 2 )(f> = 0, (5) 

which is the usual Klein-Gordon equation in Minkowski spacetime with a time- 
dependent mass term which vanishes if the field is conformal. 
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2.1. The state of the quantum field 

Let us discuss the kind of states that we will consider for the scalar field 0. The Klein- 
Gordon equation (5) admits the following complete set of orthonormal solutions: 

« £ (i7, x) = fk(v) {2n) l_ 1)/2 e iU - 3 , (6) 

with k G M. n ~ L and k = \k\. The mode functions fk satisfy the equation 

h + (k 2 + M 2 )f k = 0. (7) 

In order for the set {u^} to be orthonormal, the mode functions must also satisfy the 
Wronskian condition 

hti - fkfk = (8) 
The solution of equation (7) is unique once we specify initial conditions at some initial 
time r/i, fk(i]i) an d fkijl'i)i which must also be consistent with (8). Given a set of 
orthonormal modes, the field operator can be expanded in terms of the associated 
creation and annihilation operators as = + an d a Fock space based on 

the vacuum defined by Sg|0) = can be constructed. Therefore, since different choices of 
initial conditions for f^ give rise to different sets of mode functions, they define different 
(homogeneous and isotropic) vacua. In this paper we will proceed as follows. Define an 
auxiliary scale factor a^,(r]), with domain (— oo,^], such that 

lim Ml = 0. (9) 

1]^ — CO 

Then the initial conditions for the mode functions will be chosen as 
/*(»») = fk(Vi) 

fk(Vi) = fk(Vi), (10) 
where ffiji) is the solution of the mode equation (7) with M|(?7) replacing M 2 (rj), 

fZ + (k 2 + Ml)f* = 0, (11) 
that behaves as a standard plane wave at past infinity, i.e. 

f!M - ^ (12) 



when i] — > — oo. The state that we will consider is the vacuum associated with 
the initial conditions (10). By taking different auxiliary scale factors a^, one can get 
different initial conditions, and thus different vacua. So far, a* is completely arbitrary, 
except for condition (9). However, in order for the state to have the right ultraviolet 
behavior, further conditions must be imposed on a^,, as we will see below. 

The states that we have defined have a remarkable computational advantage: they 
are just the m-vacuum when the scale factor is 

a(rj) = { a{v ) . 1° 1 71 > Vi (13) 
j a* (77) for i] < r\\. 

Therefore, instead of working with the scale factor a(i]) and the state one can 
perform the computations regarding the scalar field as if it were in a RW spacetime 
with scale factor a{rj) and its state were the m-vacuum evolving from r\ = —00. 
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2.2. The expectation value of the stress-energy tensor 



a n ~ l J d n - l xT£. (14) 



The trace of the stress tensor operator is obtained by functional differentiation of the 
classical action (4) with respect to the scale factor: 

5S_ 

5a 

Similarly, it is well-known [52, 53] that its expectation value in the state 
can be obtained by functionally differentiating the Feynman- Vernon influence action 
S'jp [a + ,a~], which is defined as 

e iS*[a + ,a-] = J p0+D0-p v ,[0+(r ?i ),0-(r7 i )]e< 5[a+ ^ +] - 5 ^'^ ] ), (15) 

where is the density matrix corresponding to the state and the field configura- 
tions <p + and <p~ are supposed to coincide at some final time r]?. Indeed, functionally 
differentiating this influence action yields the expectation value of the trace of the stress 
tensor in the state |\&), 



5S* 



5a 



= a 



n-l 



V(T£)*, (16) 



e iSi F {a+,a ] 



where V is the volume of space, V = J d n ~ 1 x, and we have taken into account that 
the state |^) is homogeneous. Following the discussion around equation (13), instead 
of S^p[a + ,a~] we can work with Sif[o, + ,d~], where the absence of the superscript \1/ 
indicates that the state under consideration is the m-vacuum. The latter influence 
action has the simpler form 

J P0 + P0V( s[s+ '* +] - s ' [a ~ •*"]). (17) 

The complex conjugate of the classical action appears because the usual — ie prescription 
is used. Since the action (4) is quadratic in the field, this path integral is Gaussian, and 
a formal expression for the influence action is readily obtained: 

S IF [a + ,a-] = -^trmG, (18) 

where G is the inverse of the matrix A = diag(A + , A_), with A + = i]^ d^d u — M++ie, and 
A_ = — {rf v d ll d v — M 2 — ie) . The function M 2 {rf) has been defined below equation (4). 
The bar and the subscript ± just indicate that it corresponds to the scale factor a ± . 
Functional differentiation gives 

<>- : vrr;;). (19) 



5a 

where (T^) is the expectation value, in the in-vacuum, of the trace of the stress tensor 
with the scale factor d{rj). For 77 > 771 , we have 

(f£(v)) = (T£(v)W (20) 

Since the state is homogeneous and isotropic, the only independent components 
of the expectation value of the stress tensor are the trace and the 00 component (the 
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energy density). § In fact, the latter can be obtained from the former. Indeed, using 
Gauss's theorem together with the equality V M ((T V1I ')£ I/ ) = (a/a)(Tff) for £ = d/dr], 
which is a consequence of the stress tensor conservation law and the fact that £ is a 
conformal Killing vector field of the RW spacetime, we get 

a n - 2 (v)(Toom = ~ f W^\rf)t(Tf)(Tj!(Tf)), (21) 

J — OO 

where we took into account that for a field in the in- vacuum d n ~ 2 (77) (T 00 (r))) vanishes 
in the limit i] — > — oo. Again, for rj > rji we have 

(foofa)) = (Too(»?)>». (22) 
2.3. Renormalization 

The influence action (18) suffers from ultraviolet divergences in the limit n — > 4. The 
divergent contributions to the influence action can be written as 5di V [g + ] — Sdi v [g~], 
where <Sdi V [<7] corresponds to the divergent part of the in-out effective action. The short- 
distance behavior of the in-out effective action of a free scalar field has been thoroughly 
studied for a general metric g (not necessarily of RW type) and it has been shown that 
the divergent part has the following form [40]: 



n — 4 

S *' |S| = - 324,-4) ' d " X ^ 



m 4 + 2m 2 vR + v 2 R 2 



+ ^ (R^R^ - R^Br) ] , (23) 

where /i is an arbitrary mass scale introduced for dimensional consistency which plays 
the role of the renormalization scale when using dimensional regularization; R^ ua /3 and 
Rfj, u are the Riemann and the Ricci tensors respectively. The dynamics of the metric in 
semiclassical gravity can be derived from the following closed-time-path (CTP) effective 
action [52]: 

r CT p[s + , Si = S g [g + ] - S g [g-\ + S 1F [g + ,g-\, (24) 

where S g [g] is the classical gravitational action, which will be fully specified below. 
Physical predictions can still be finite provided that the so-called bare parameters of 
the gravitational action have the appropriate divergent behavior so that they cancel 
the divergences of the bare influence action Sip[g + , g~] and the total effective action 
^CTp[g + , g~] is finite. This is accomplished by considering a bare gravitational action 
of the form 

S g [g] = Sl en [g] - S div [g], (25) 

where the finite renormalized gravitational action is given by S^ 011 ^] = ^ehM + S c [g]. 

The first term is the Einstein-Hilbert action 

1 



SehM = J d'x^g- (R - 2 A) , (26) 

§ Strictly speaking the energy density physically measured by a comoving observer corresponds to T tt , 
which is related to Too by T tt — a~ 2 Too- However, from now on we will loosely refer to Too as energy 
density as well. 
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where A is the cosmological constant and lp is the Planck length, lp = y/G, with G 
being the gravitational constant. On the other hand, S c [g] corresponds to the additional 
counterterms, which are quadratic in the curvature, and is given by 



S c [g] = J d A x^~g [a (R^R^ - R^) + PR 2 ] , 



(27) 



where a and j3 are dimensionless parameters. Note that since the bare parameters should 
be independent of the renormalization scale fi, the renormalized parameters 1/G, A/G, 
a and (3 depend on ix in such a way that they cancel the \i dependence of the finite terms 
that arise when taking the n — > 4 limit in (23). More specifically, under a change the 
renormalization scale //—>// each one of the renormalized parameters changes by a term 
proportional to ln(/////). This also implies that the total effective action TcTp[g + , g~] is 
invariant under the RG (i.e. independent of /i). Taking (25) into account, the effective 
action in (24) can be written as r CTP [g + , g~] = S™ n [g + ] — S™ n [g~] + Slp n [g + , g~], 
which is manifestly finite and where S{p n [g + , g~] = Sif[<? + , g~] — Sdi v [g + ] + £div[<7 - ]- 
Moreover, to make some of our expressions below more compact, we will reabsorb the 
part of the renormalized gravitational action which is quadratic in the curvature into 
the renormalized influence action so that we have 



sing" 



Sdw[g + ] + S div [g ] + S c [g + ] - S c [g~ 



(28) 



,9 ] = Si F [g + ,g~ 

and the total effective action becomes 

r C Tp[s + ,s1 = S E u[g + ] - S EH [g-} + S^\g + ,g-]. 

Let us specialize these general results to the particular case that we are considering 
For a metric of the form (1), one has 

2 / . \ 4 



(29) 



/ 



d n xV=g (R^ P aR^ pa - R^Rn = - V(n - 4) J d V 

+ 0((n-4) 2 ). 



(30) 



where we neglected a total divergence that does not contribute to the equation of motion. 



Using (3) and the fact that ^f^g 



[l + (n- 4) In a) + O ((n - 4) 2 ) one can 



proceed analogously for the terms in the integrand of (23) which are proportional to 
R 2 , R and a constant. Substituting into (28) and taking into account the definition 
of the time-dependent mass below equation (4), we obtain the following result for the 
renormalized influence action of our scalar field model: 

... 



o 



IF 



a , a 



S 



if a 



+ 



32vr 2 



A 



dr]{ ( - 



+ M 4 



1 
90 



1 



n - 4 



+ ln(a/x) 



(31) 



where we have introduced the difference notation, Af = f + — f~ , and some numerical 
factors have been absorbed into the free parameter /3. The 0(n — 4) terms have been 
neglected, as they vanish in the limit n — * 4. Note that since the right-hand side 
of (30) vanishes when n — > 4, our renormalized influence action is independent of 
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the renormalized parameter a. On the other hand, the logarithmic dependence on 
/i is compensated by the \x dependence of (3 as well as 1/G and A/G in the Einstein- 
Hilbert action so that the total effective action is independent of //, as explained earlier. 
Moreover, although the argument of the logarithm in (31) has dimensions of mass, 
everything becomes well-defined when this term is combined with the bare influence 
action, as we will see below. From this renormalized influence action, one can obtain a 
finite expectation value of the stress-energy tensor by following the steps described in 
the previous subsection. 

We close this subsection with a few remarks about the renormalized gravitational 
coupling constant G ([/,), the renormalization scale \i and the experimental value of G 
measured with a Cavendish-type torsion balance. In this sense, an important question 
that needs to be addressed is how the result of this kind of experiment can be used to fix 
the value of G(fi) at some scale \x and what is the most natural scale to consider. The 
Cavendish experiment can be essentially understood as a two-body scattering problem 
in the infrared limit. This is possible provided that there is a large separation of scales 
between the typical radius of curvature of the cosmological background (and in practice 
other local gravitational fields) and the Compton wavelength of the lightest massive 
particle: we need the typical scale of the experiment to be much larger than that 
Compton wavelength while being able to subtract or neglect other contributions to 
the local gravitational field. The condition on the Compton wavelength of the lightest 
massive particle guarantees that we are in the decoupling limit for all the massive 
fields, whereas massless fields do not contribute to the renormalization of 1/G and their 
correction to the Newtonian gravitational potential is suppressed by a factor (Zp/r) 2 , 
with r being the distance between the two masses (the typical scale of the experiment) 
[54, 33]. In that case, one can see that for a renormalization scale of the order of 
the mass of the lightest massive particle, /i ~ m, G(/i) corresponds to the measured 
value if there are no additional massive fields. Otherwise, the measured value of 1/G 
would correspond to 1/G(m) plus a contribution of order m\\n{m / m^) for each heavier 
field with mass m^. Note that unless the mass m of a field is not too far from the 
Planck mass, its contribution to the running of 1/G(fi), which is of order m 2 ln(/x/ //) is 
extremely small compared to the value of 1/G(fi) itself. On the other hand, the value of 
the cosmological constant, which is much more sensitive to the running due to massive 
fields (even light ones), should also be determined in a similar way from cosmological 
observations. Finally, it is much harder to place upper bounds on the values of a and (5 
because the corrections from the terms quadratic in the curvature are very suppressed 
for phenomena with typical scales much larger than the Planck length. 

2.4- The semiclassical Einstein equations 

The semiclassical Einstein equation can be derived from the CTP effective action, 
as given by (29), by functionally differentiating with respect to g + and then taking 
g + = g~ = g. The Einstein-Hilbert action gives the Einstein tensor plus a cosmological 
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constant term, and the renormalized influence action gives the renormalized expectation 
value of the stress tensor operator [52] : 

G^ + Ag lu , = 8'Kll{T lu ,)v, (32) 

where G^ v is the Einstein tensor. Once the expectation value of the stress tensor is 
known, the backreaction of the quantum field on the mean geometry can be computed by 
substituting it into the semiclassical Einstein equation. Because of spatial homogeneity 
and isotropy, only the equations for two components are independent. They can be 
chosen to be, for instance, the 00 component and the stress tensor conservation law. 
Since the latter has already been taken into account when deriving the energy density 
from the trace, we only need the 00 component, which is the semiclassical analog of the 
Friedmann equation, and in four dimensions has the form 

a 2 = FV + ^a 2 (r 00 ),, (33) 

where we have introduced the Hubble constant H = a/A/3. Note that, since we only 
know (T 00 )^ for r\ > r]i, this equation can only be solved for r] > r^. Two time-scales 
appear in the equation: the Hubble time H' 1 and the Planck time tp — lp. In order 
for the semiclassical approximation to be valid, these two time-scales have to be well 
separated, i.e. we need l P H <C 1. 

As we will see below, the expectation value of the energy density contains up to 
third order time derivatives of the scale factor. Therefore, equation (33) is a third order 
(integro-)differential equation, and its space of solutions is much larger than the one 
corresponding to the classical Friedmann equation. However, many of these solutions 
have characteristic time-scales of the order of the Planck length and lie beyond the 
domain of validity of our low-energy EFT approach. In the spirit of this EFT approach, 
we will look for perturbative solutions in powers of (Hip) 2 : 

a( V ) = a Q ( V ) + {l p HfaM + O ((/ P #) 4 ) . (34) 

A perturbative expansion may sometimes miss the right long-time behavior of the 
semiclassical solution. This can happen when the effect of the quantum corrections, 
although locally small, builds up over long times giving rise to substantial deviations 
from the classical solution. An example of such a situation is the evolution of a 
black hole spacetime when the back reaction of the emitted Hawking radiation is 
taken into account. One possibility in those cases is to modify the backreaction 
equation using an order-reduction procedure and then solve the resulting equation 
non-perturbatively [32, 55]. If these cumulative effects were present in our case, the 
correction a\ would get large and the perturbative approximation (34) would no longer 
be valid. However, our perturbative treatment would still be useful to signal that these 
cumulative effects are taking place. In fact, one can check the validity of the perturbative 
approximation a posteriori by making sure that the perturbation remains small at all 
times. 
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Substituting (34) into the Friedmann equation (33) and solving order by order in 
perturbation theory, we get to lowest order that a is just the scale factor of de Sitter 
spacetime, 

aoiv) = adsiv) = -j^j, (35) 
with i] < 0. The first order equation is 

ai = -\ + ^(T 00 )f, (36) 

where (T 00 )^, s is the expectation value of the energy density evaluated on the de Sitter 
background. This is a first-order linear differential equation whose solution is readily 
obtained: 

aM = 3^f f d n'v ,2 (TM))f- (37) 

In equations (35) and (37) we have dropped an arbitrary integration constant, because 
it can be set to zero by an appropriate choice of the origin of conformal time (see [50] 
for details). 

In the next section we will compute the expectation value of the energy density 
for non-conformal fields in a general RW background according to the steps described 
in this section, and then we will specialize it to a de Sitter background to obtain the 
first-order semiclassical correction to the de Sitter spacetime from equation (37). 

3. Weakly non-conformal fields 

In this section we will obtain an explicit result for the influence action Sip[d + ,a~], 
making use of the formal expression (18) and the assumption that the field is weakly 
non-conformal, that is, M 2 /a 2 -C H 2 . The expectation value of the stress tensor in the 
state will then be derived and the semiclassical equation solved for the particular 
cases (c = 0,m^0) and (v ^ 0, m — 0). 

For a weakly non-conformal field, the matrix G in equation (18) can be computed 
perturbatively in powers of M 2 . Following [56, 57, 58] we define A = A + V, where 
the matrix V includes the time-dependent interaction with V + = —M^, and V- = M 2 . 
Then, up to second order in M 2 , 

G = G°(1-VG° + VG°VG° + ...), (38) 

where G° = (A )' 1 is the Minkowski, massless CTP propagator. Its components are 
G° ++ = A F , G°„_ = -A D , G° + _ = -A+ and G\ = A", where A F and A D are 
respectively the Feynman and Dyson propagators and A ± are the Wightman functions: 

f d n p e ip ' x 

A V°W=-J (27T)V^' 

A ± (x) = ±2ni J ^e^5(p 2 )6(TP°). (39) 
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Substituting into equation (18) we have (up to second order in M 2 ) 
S lF [a + ,a~} = - ^trlnG - ^tr(M 2 A F ) - ^tr(M!A D ) 

- l -ti{M 2 + A F M 2 + A F ) - l -ti{MlA D M 2 _A D ) 

- -tr(M 2 A+M 2 A~). (40) 
2 

The first three terms in this equation do not contribute to the expectation value of the 
stress tensor. The first one is independent of a, whereas the second and third terms are 
tadpoles which are identically zero in dimensional regularization [59]. Therefore, the 
influence action is quadratic in M 2 . Computing the last three terms and expanding in 
powers of (n — 4), we get 

Re * IF = -3^/^U 



-^AM 4 (t?) + 2AM 2 (t?) J dr)'H(r) - rf; #)£M 2 (t/) j , (41) 

where the terms of order 0(n — 4) have been neglected. The range of both 77 and 77' 
integrations is (—00,77/], where rjf is an arbitrary final time (larger than any other time 
that we may be interested in). We have used the difference and semisum notations, 
A/ = /+ - /- and £/ = (/+ + f~)/2 respectively. As shown in [57], the kernel 
H(rj — 77'; 9) is given by 

H(v ~ 1/; 9) = J ^ + f sign(-c)) , (42) 

where 9 2 = exp(2 + ln47r — 7), and 7 is the Euler-Mascheroni constant. In equation (41) 
we have dropped an imaginary term, that does not contribute to the expectation value 
of the stress tensor because it is proportional to (AM 2 ) 2 . This is consistent with the 
fact that the expectation value of an observable must be a real number. Using equation 
(31) and taking the limit n — > 4, we obtain the renormalized influence action: 



2 , / x \ 4 



+ 



A(M\r]) In afa)) - 2 AM 2 (77) J drfH(rj - r]'; n)Y,M 2 {r]')\ , (43) 

where the numerical factor 9 has been absorbed into /_t. Notice that as has been 
anticipated in subsection 2.3, the term involving ln/i in (31) combines with the term 
with In \oo\ in (42), rendering the argument of the logarithm dimensionless. 

Let us finally turn our attention to the non-local part of the influence action. The 
Fourier transform (42) can be explicitly computed. If we define the functional 

K\f; V )= / drfH( V -rf^)f(rf), (44) 



using equation (VII. 7. 18) in [60], it can be shown to be 



K[f;r ] ) = -\im \r ' -^/(i/) + (lne + ln/i + 7)7(77) j> . (15) 
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In order for the integral in this equation to be finite in its lower limit, / must vanish 
at past infinity. Assuming that this condition is satisfied, the functional «[/; 77) is finite 
at i] provided that / is continuous at rj. Indeed, in this case near rj' — r] we can write 
/(?/) — f(v) + 0(r) — rj'), so that the integral in (45) has the form 
r r '~ e dn' 

/ —fW) = -f(v) ln e + finite terms, (46) 

J-oc V~V 

and the contributions that diverge in the limit e — > + cancel. Therefore, in order for 
the functional «[/; 77) to be well-defined, the function / must be continuous and vanish 
when rj — > —00. On the other hand, for a function satisfying these conditions it is easy 
to see that 

±K[f;ri) = K[f;v)- (47) 

The expectation value of the trace of the stress tensor, which will be obtained by 
functional differentiation of the influence action (43), will include a term proportional to 
[M 2 ;?]), and in some cases there will also be terms proportional to its first and second 
derivatives, as will be shown below. The condition (9) already ensures that M 2 vanishes 
at past infinity. Now we see that M 2 must also satisfy some "smoothness" conditions. 
In particular, in order for the trace not to diverge at the initial time rji, the scale factors 
and a in equation (13) must be matched smoothly enough. This corresponds to 
choosing regular initial states, namely, adiabatic vacua of the appropriate order, as is 
explained in the appendix. 

3.1. Massive conformally coupled fields 

Here we restrict our attention to the particular case of conformally coupled fields (with 
v — 0), but keeping m/0. In this case, the condition (9) implies that vanishes at 
past infinity. The renormalized influence action (43) takes the form 

Re ^ IF= 3^/*^ (!) + 95 A GD + m4 [ A (« 4ln «)-2Aa 2 K pa 2 ]]}, (48) 

where we have dropped the superscript "ren" for simplicity, and k is the functional 
defined in equations (44) and (45). Functionally differentiating with respect to a + and 
making use of equation (19), we obtain the expectation value of the trace of the stress 
tensor, which after some manipulations reads 



1 • /^,m d [ T . d ( a 



1 / a 
2 



960tt 2 



-\ 4 

a 



4 ,_ 

m 



;l <r\na ■ l'<r_y (49) 

where we again absorbed some numerical factors into the free parameter f3, and the 
functional F is defined as 

P dn 2 

F[a 2 ;r ] ) = -2j drf^n^r/). (50) 
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We see that the expectation value of the trace of the stress tensor involves a term 
proportional to K,[a 2 ]. According to the discussion below (45), in order for the trace to 
be finite at the initial time rji, the auxiliary scale factor employed in the construction 
of the vacuum must verify 

a*(Vi) = (51) 

Using equations (21) and (22) it is straightforward to obtain the expectation value of 
the energy density for rj > rji given a scale factor airf) and an initial state |\&): 



a 2 {T 00 )y = (3 



d 



drj \a 2 



1 (a 

2\a 



+ 



1 



960tt 2 



/ a 
V a 



m 



32n 2 



a \na + F[a 



(52) 



This result is valid for any scale factor a(i]) and for any vacuum \^/) (as defined in 
subsection 2.1) satisfying the condition (51). Let us now particularize this to the case of 
a de Sitter background. Since de Sitter spacetime has constant curvature, the first term 
in (52) vanishes: its origin is the term proportional to R 2 in the influence action, which 
yields a term proportional to OR after functional differentiation. Therefore, there is no 
dependence on the free parameter (3, and one is left with 



ds (^oo)* S 



m 



960tt 2 (-r]Y 32tt 2 



-H V y 



ln(-#7/) + F[477) 



where a 2 , (77) is defined as 



'dS 



for i] > rji 
for i] < r/i. 



a\ — a^g + b~(i\ 

Note that, as a consequence of the conditions (9) and (51), we must have 



lim Sa^, 



0, <$4(t&) = o. 



(53) 



(54) 



(55) 



One possibility is, of course, 5a%{ri) =0 M rj < r/i] in this case, the state is the m- vacuum 
of de Sitter spacetime, which is the so-called Bunch-Davies vacuum. Substituting the 
scale factor (54) into equation (50), we get after some manipulations that the state- 
dependent part of the expectation value of the energy density can be written as 

F[^-r,) = F[a 2 dS -,r,) + F[8al-,r h ) 

drj" 



-2 dr/5a%(rf)V f 4z 

J —CO J —00 V 



^ a dS 



rj" dr] h 



(56) 



where V denotes Cauchy's principal value. The first term in this equation, as well as 
the second integral in the third term, can be explicitly computed. The latter involves a 
logarithm that can be expanded in a Taylor series. The result is 



F[al;rj) 



-HrjY 



In(-iir)) + 7 - 



-Y. A n(- H v) n - 2 -B*M-H V ), 



(57) 



n=0 
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where the state-dependent, dimensionless constants A% and By are given by 

^7F(^/_1(=^4(,') for -^2 
4 r Vi dn' 

^/If^ 4 ^ (58) 

All these integrals are finite due to the first of the conditions (55). Note also that they 
all vanish when the state is the Bunch-Davies vacuum. Substituting the result (57) into 
equation (53) we obtain an explicit result for the expectation value of the energy density 
of the scalar field in de Sitter spacetime for the state |^): 

'Y J At{-H V ) n + B^-H V f\n{-H V ) , (59) 



(7oo)* S — (^oo)bd + 



327T 2 



_n=0 

\dS 



where (T o) B f) is the same expectation value in the Bunch-Davies vacuum, 



(Too) 



dS 
BD 



{-HtjY 



H m f fi 



In — + 7 - 



(60) 



_960tt 2 32tt 2 V H ' 4 / 

This result agrees through order m 4 with the exact result in equation (6.183) of [40]. 
To show the equivalence one needs to take into account that their renormalization 
prescription corresponds to taking fi equal to m times a dimensionless number involving 
the Euler-Mascheroni constant, and the fact that their counterterms are given by (23) 
but with (f - 1/6) instead of v = (f - Q = (f - 1/6) - (n - 4)/36 + 0((n - 4) 2 ). Note 
also that the /x-dependent term in (60) is exactly compensated by the \i dependence of 
the cosmological constant in the renormalized Einstein-Hilbert action. || Moreover, there 
is no /i dependence in the state-dependent part of (59), as can be seen from examination 
of the coefficients (58). The only one that could potentially depend on fi is A* , but 
it actually does not: the /z-dependent part of K,[f; if) is local, as can be seen from its 
definition (45), so that the /x-dependent part of F[5a\;rjj) is proportional to 



/Vi 
dr]'5a% 
-oo 



d{5a\. 



= (Km) = o, 



dr]' 

as follows from conditions (55). Therefore, the semiclassical Einstein equation is ji- 
independent, as emphasized in subsection 2.3. 

Another point that one must check is that the infinite series in (59) converges. To 
see that this is indeed the case, note that the absolute value of the coefficients A% in 
(58) (for n > 2) is bounded from above as follows: 

^l^ 4m ^l^l' n(n-2)W - <61) 

|| There is no u dependence associated with the term proportional to the Ricci scalar because there is 
no divergence associated with that term for v = 0, as can be seen from (23). However, for the massive 
and non-conformally coupled case there would be such a dependence. 
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Therefore, the absolute value of the n-th term of the series is 



|A*|(-#77)"<4max(|54| 



1 



< 4 max(|<5a|| 



n(n — 2) 



,(62) 



n(n - 2) 

where the last inequality holds for rj > w\ (recall that both rj and 771 are negative). Now, 
the infinite series in (59) converges because the series ^ 1/n 2 , which is Riemann's zeta 
function C(2), is finite: 

00 1 2 



As a final remark, notice that the expectation value (59) tends to the Bunch-Davies one 
(60) at future infinity, 77 — > 0, regardless of the state 

The final step is to use the result (59) to compute the semiclassical correction to 
the scale factor of de Sitter spacetime using equation (37). The result is 

C*(-H V r + Dq,(-H V ) 3 \n(-H V ) 



01(77) 



(v)- 



m 



2AnH 4 



(64) 



where af D (7?) is the correction associated with the Bunch-Davies vacuum, 



.bd 



m 



7207T 247T# 4 



(65) 



and the state-dependent coefficients and Dq, are 
A? 



n + 2 



for n 7^ 3 



4* 

n - B * 

— ——. 

5 



25" 



(66) 



Substituting this result into the expansion (34), we finally conclude that the solution of 
the semiclassical Friedmann equation (33) is 



a( V ) = - -I- - (l P Hf 
Hi] 



m 



n=l 



24irH A 

+ O ((/ P iJ) 4 ) , 

where the corrected Hubble constant H is given by 

1 



J2 C n(-Hv) n + D^-HrjfH-Hrj) 

(67) 



H — H <1 + (IpHY 



m 



1 A 4 3 

In — + 7 

H 4 



(68) 



_720vr 2A-kH a 

Note that the explicitly /i-dependent term in (68) is compensated by the \x dependence 
of H, which it inherits form the \x dependence of A, so that H is /i-independent. [The 
/j, dependence of H in the terms of order (IpH) 2 is irrelevant since it would give rise to 
contributions of higher order in (IpH) 2 .] If the state chosen is the Bunch-Davies vacuum, 
then the coefficients and Dq, vanish, and the solution (67) of the semiclassical 
Friedmann equation is de Sitter spacetime with the corrected Hubble constant H. If 
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we choose any other state, then the solution is no longer the scale factor of de Sitter 
spacetime, but tends to it at late times (as rj — > 0~). Finally, notice that since there 
is no dependence on the free renormalization parameters a and /3, this result is fully 
predictive once the cosmological constant and the gravitational coupling constant (or, 
equivalently, /p) are measured. 

It should be emphasized that the logarithmic term which appears explicitly in (53) 
and comes from the local term proportional to m 4 and In a in (48) is cancelled out by 
the first term on the right-hand side of (57), which contains the contributions from the 
non-local term in (48). If it had not been cancelled out by the contribution from the 
non-local term, such a local term would have dominated the result for the stress tensor 
at late times, precluding the existence of a self-consistent de Sitter solution and altering 
our main conclusions about the semiclassical stability of de Sitter spacetime. This is an 
important point whose implications will be further discussed in section 5. 



3.2. Massless non-conformally coupled fields 



Next we move to the opposite situation: we will consider massless fields but non- 
conformally coupled (i.e. m = but v ^ 0). We will proceed in strict analogy with the 
previous subsection, so some of the details will be skipped (this case has been considered 
in great detail in [50]). In this case, M(rj) in four dimensions reads 



M 



6z/— 

a 



Quu, (69) 

as follows from equation (3). Condition (9) translates here into uq, — > at past infinity. 
Specializing the influence action (43) to this case, functionally differentiating it with 
respect to a + and applying equation (19), one finds that the expectation value of the 
trace of the stress tensor depends on the functional k[u; rj), as well as its first and second 
derivatives. Therefore, from the discussion around equations (46) and (47), we see that, 
in order to keep the trace finite at the initial time, the scale factor a* employed to 
generate the state \$>) must satisfy the following conditions: 

uvim) = u (vi) = u(m) = u(m)- ( 70 ) 

In other words, the scale factors before and after the initial time must be matched with 
continuity up to the fourth derivative. Once the trace is known, the expectation value 
of the energy density of the field for rj > rji, given an initial state and a scale factor 
a(r/), is obtained from equations (21) and (22). The result is 

2" 



a 2 (Too)* = P 




+ 



1 



960tt 2 



d 2 a 



where the non-local, state-dependent part reads 



— + T[a,u}\, (71) 



dr)'u(r)')K[ii;r)'). (72) 
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Note that the first and second terms in equation (71) are just the same as in (52). They 
both come from the first two terms in the influence action (43), which are independent 
of the particular form of the time-dependent mass M(i]). In order to specialize this 
result to a de Sitter background, but keeping the state general, one must compute 
T[a dS ,u ;v) witn 

- / x J MdS for 7] > T)i 

! Uq, = U dS + OUq, for 1] < 

where u d s(v) = 2/r^ 2 [u^(i]) was denoted by v{rf) in [50]]. Note that conditions (9) and 
(70) translate into the following conditions on Su^: 

lim Siiif, = Su^^rji) = 5u^(r]i) = 5u^(r]i) = 0. (74) 

Once again, taking 5uq, = corresponds to the in-vacuum for de Sitter spacetime, i.e. 
the Bunch-Davies vacuum. Substituting the scale factor a d $ and the function uq into 
equation (72), we get that the state-dependent part of the expectation value of the 
energy density in de Sitter spacetime is 



T[a dS ,u ;r}) = -H 4 



n=0 

where the dimensionless state-dependent constants P* and are given by 

n 2 + n-2 [ Vi dr]' 



,(75) 



4 5 

In(-Hr)') - + k[6uv, rf) 



-if y " {-if 



Q ^^Lrw SuM) - (76) 

These state-dependent coefficients are finite, as a consequence of the first condition in 
(74), and vanish in the particular case of the Bunch-Davies vacuum. Making use of this 
result, it is now straightforward to explicitly obtain the expectation value of the energy 
density for the state |^) in de Sitter spacetime: 

9z/ 2 P 4 



(Poo)* — (Poo)bd + 



4vr 2 



Pn(-Hv) n + Q^-Hvf H-HV) 



,(77) 



_n=0 

where the first term is the expectation value in the Bunch-Davies vacuum, and reads 

(Too)b S d = 7c|^ (gg^ " ^) ' (78) 

In this case, there is no /i dependence at all. First, there is no dependence in the 
state-dependent part: the only coefficient which could depend on /j, is P 2 * and one 
can show that it actually does not, exactly in the same way as done for A* in the 
previous subsection. Second, there is no dependence on \x in the state-independent part 
either because the only term that depends on fi is proportional to OR, which vanishes 
when evaluated on a de Sitter background. This term is related to the PR 2 term of 
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the renormalized gravitational action. On the other hand, there is no dependence on ji 
associated with the terms of the Einstein-Hilbert action because there are no divergences 
involving those terms in the massless case, as can be seen from (23). 

The convergence of the infinite series in (77) can be checked by following the same 
steps as in the previous subsection. One also needs to use the relation 

drf , , A 1 [ Vi H 



J 

J — c 



[-r,')n+i 



5uy(r)') = 1 f dT } 8un,{ri), 
n(n - 1) J_ 00 (-?/)" 1 



which is a consequence of the conditions in (74). Another interesting point is that, 
analogously to the massive conformally coupled case, the expectation value (77) tends 
to the Bunch-Davies one in the asymptotic future (i.e. when rj — > 0~) regardless of the 
state |^). 

We can now compute the semiclassical correction to de Sitter spacetime, making 
use of the expectation value (77) and equation (37). Substituting the result into the 
expansion (34), we obtain that the solution to the semiclassical Friedmann equation is 



yZ R n(- H v) n + S^(-HrjfH-Hrj) 



n=l 



Hi] vr 

+ O ((/ P # ) 4 ) , (79) 

where the relation between the coefficients in the stress tensor and the solution for the 
scale factor is the same as in the previous subsection: 

Rl = S± for n ^ 3 
n n + 2 r 



P-2 _ 



S* = %. (80) 



In this case the corrected Hubble constant is 



H = H 



(81) 



720% 2tt J 

If the state chosen is the Bunch-Davies vacuum, the solution to the semiclassical 
Friedmann equation is just de Sitter spacetime. For any other state, the solution is 
the scale factor of de Sitter spacetime plus some terms, but these terms vanish at 
future infinity. The result (79) is fully predictive, as there is no dependence on the free 
renormalization parameters a and (3. 



4. Strongly non-conformal fields 

In this section we briefly describe how to obtain the influence action for strongly non- 
conformal fields corresponding to a large mass m 2 ^> H 2 (the curvature coupling 
parameter v can take arbitrary values as long as m 2 ^> vR). 

Whenever the mass of the fields is much larger than the inverse of the curvature 
radius, one can introduce a quasi-local expansion for the Feynman propagator and 
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other related Green functions (see section 3.6 of [40]). This can then be used to 
calculate the influence action as a local expansion of terms involving positive powers of 
the curvature times the appropriate power of the mass squared which gives the right 
dimension, as explained in sections 6.1 and 6.2 of [40]. This kind of expansion is known 
as the Schwinger-DeWitt expansion and it is often employed within some covariant 
regularization scheme in order to identify the divergent contributions to the effective 
action, which can be dealt with by a suitable renormalization of local terms in the bare 
gravitational action. When identifying divergent terms in four-dimensional spacetimes, 
the expansion is truncated beyond order one and only includes terms of order m 4 , m 2 and 
1. After subtraction of the divergences associated with those terms through the usual 
renormalization procedure, one is finally left with the same kind of terms multiplied by 
finite coupling constants. There are, however, additional finite contributions from terms 
in the expansion corresponding to negative powers of m 2 . In particular, the terms of 
order m~ 2 were explicitly considered in [61], where the resulting effective action was 
finally used to calculate the expectation value of the stress-tensor operator for fields in 
a Kerr black hole spacetime. 

In fact, using arguments based on the relevant symmetries (diffeomorphism 
invariance in this case) and power counting, one can easily infer the form of the influence 
action in this regime except for the precise values of dimensionless coefficients. This 
kind of approach is very common when dealing with effective theories and it is nicely 
illustrated by the example of the Euler-Heisenberg Lagrangian [62], which describes the 
effects of electron loops on photons when the energies of the photons are much smaller 
than the electron mass. Since the influence action is described in the regime m 2 ^> R by 
a local expansion, it can be separated into two parts: Sfp n [g + ,g~] = S c s[g + ] — S e a[g~]. 
S e s[g] coincides with the result of integrating out the massive fields in the in-out 
formalism and it is given by 



where a, b, c and d are dimensionless parameters, C iiVpu is the Weyl tensor and 0(R 3 /m 2 ) 
denotes terms of order 1/m 2 or higher, which involve cubic or higher powers of the 
curvature but also terms with covariant derivatives starting with those quadratic in 
the curvature and with two covariant derivatives. Note that terms corresponding to 
a total divergence (such as the DR term) or a topological invariant have not been 
included. That is also the reason why only two of the three possible terms quadratic in 
the curvature have been considered (the third one being R^R^) since there is a linear 
combination of the three which corresponds to the Gauss-Bonnet topological invariant. 
The first four terms that appear on the right-hand side of (82) are in general multiplied 
by divergent coefficients in the bare influence action; % the renormalized influence action 
is rendered finite by including suitable counterterms in the bare gravitational action, 
as explained in section 2.3. (The parameters a and b exhibit a fi dependence which is 

% However, the divergent coefficients of the Ricci scalar and a certain linear combination of the R 2 and 
(jtivpa 'C terms vanish when v = 0. 




<Pxy/=g [am 4 + bm 2 R + cR 2 + dC^C^ + 0(R 3 /m 2 



)]. 



(82) 
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compensated by the \x dependence of 1/G and A/G in the renormalized Einstein-Hilbert 
action.) There is, nevertheless, an inherent ambiguity in the finite part of the coupling 
constants for those counterterms, which can only be determined experimentally. On 
the other hand, that is not the case for terms of order 1/m 2 and higher, which do not 
need to be renormalized and can be explicitly computed using a quasi-local expansion 
as described above. The value of those coefficients can then be determined exactly 
assuming a vanishing contribution from the bare gravitational action, and even if there 
are contributions from physics at higher energy scales, they will be highly suppressed 
provided that there is a sufficiently large separation of scales. 

By functionally differentiating (82) with respect to the metric one can obtain 
the expectation value of the stress tensor operator for the massive field, which is 
clearly conserved since SW^] is diffeomorphism invariant. The first two terms simply 
correspond to a finite renormalization of the gravitational coupling constant and 
the cosmological constant, which is already reabsorbed in their measured values. 
Furthermore, the terms quadratic in the curvature give no contribution when considering 
the backreaction of the massive field on the dynamics of small perturbations around de 
Sitter spacetime. This is because when using the perturbative treatment described 
in section (2.4), one needs to evaluate the stress tensor on the classical de Sitter 
background, and the functional derivatives of terms quadratic in the curvature, whose 
results can be found in [40], vanish when evaluated for a de Sitter metric. Hence, through 
the order considered in (82) there is no backreaction effect due to very massive fields in de 
Sitter. The solutions of the semiclassical equation for spatially isotropic perturbations 
simply correspond to a de Sitter spacetime associated with a cosmological constant 
corrected by the finite renormalization due to the first term on the right-hand side of 
(82). If this renormalization is already included in the classical value of the cosmological 
constant, all the perturbed solutions are just time translations of the classical one and, 
thus, physically equivalent. 

Several remarks about the state of the fields are in order. Although there is 
in general no unambiguous and preferred notion of vacuum and particle in curved 
spacetime, the quasi-local approximation for a field with a mass sufficiently large 
compared to the curvature defines a distinguished family of states which correspond 
to a natural choice of vacuum. For a RW spacetime they coincide with the notion of 
adiabatic vacua of infinite order that can be introduced in that case, and the number 
of particles of one state with respect to another state of the family is exponentially 
suppressed by a factor roughly of order exp(—m/H). If one starts with one these 
states and considers its evolution in a spacetime region where the condition m ^> H is 
satisfied everywhere, it will remain within the same class of vacua and, thus, essentially 
unexcited. On the other hand, if the evolution involves a region where the condition 
m ^> H is not satisfied, the state in the region where the condition is satisfied will in 
general be excited. However, it can still be described in terms of the natural notions 
of vacuum and particle excitation associated with the adiabatic vacua of infinite order. 
The backreaction of the field on the spacetime geometry can then be understood as a 
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combination of vacuum polarization effects and the contribution from the energy density 
of the created particles. While the backreaction due to the vacuum polarization effects is 
properly taken into account by (82) and the methods described in this section, the effect 
of the created particles can be simply described in terms of a non-relativistic perfect 
fluid (at least for free fields). 

The RG running of the cosmological constant (sometimes together with an RG 
running of the gravitational coupling constant) due to massive fields in the decoupling 
regime (corresponding to m 2 ^> H 2 ) has been proposed as a natural mechanism 
leading to a time-dependent dark-energy component [63, 64]. Those studies relied on 
general arguments rather than specific models for the quantum fields generating the RG 
running. On the other hand, the method outlined in this section gives an explicit 
result for the effective action governing the dynamics of the background spacetime 
which includes the backreaction effects of a specific microscopic model for the massive 
quantum fields. As explained above, up to the order usually considered in this kind 
of studies (the corrections from higher orders are expected to be very small) the only 
effect is simply a constant renormalization of the gravitational coupling constant (and 
the cosmological constant), which is already reabsorbed in their measured values. It 
would be interesting to see what kind of microscopic model for the fields would give 
rise to the behavior considered in [63, 64]. While the explicit calculation based on 
the Schwinger-DeWitt expansion described above was done for free fields (although it 
could probably be extended to interacting ones), the effective field theory argument 
based on diffeomorphism symmetry and power counting should apply more generally to 
interacting theories as well (as far as their gravitational interaction in the decoupling 
regime is concerned). 

5. Discussion 

In this paper we have computed the one-loop vacuum polarization for non-conformal 
scalar fields in a general spatially-flat RW background. We obtained explicit analytical 
results for two different regimes corresponding to weakly and strongly non-conformal 
fields. The weakly non-conformal case corresponds to both a small mass with m 2 -C H 2 
and a small curvature-coupling parameter u, which characterizes the deviation from the 
case of conformal coupling to the curvature. In this case we introduced a perturbative 
expansion in terms of m 2 / H 2 and u, and obtained an exact expression through quadratic 
order for the influence action, from which the quantum expectation value of the stress 
tensor operator for a family of Gaussian states can be derived. This has then been 
applied to studying the evolution of spatially-isotropic perturbations around de Sitter 
spacetime when the backreaction due to the non-conformal field is self-consistently 
included, which corresponds to solving the equations of semiclassical gravity, and we 
have calculated explicitly the solutions for all times, which are given by (67) and (79). 
There is a self-consistent solution, associated with the Bunch-Davies vacuum for the 
quantum field, with an effective cosmological constant slightly shifted from its classical 
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value due to the vacuum polarization effects. Furthermore, we have found that this 
solution is stable under spatially-isotropic perturbations since the perturbed solutions 
tend to it at late times. + On the other hand, for strongly non-conformal fields with 
a large mass the influence action for an arbitrary geometry can be expressed as a 
local expansion in terms of powers of the curvature over the mass squared, and it is 
independent of the initial state of the field. The lowest order terms simply correspond 
to a constant renormalization of the cosmological constant and the gravitational coupling 
constant which is already reabsorbed in their observed values, whereas the contributions 
from the higher-order terms only give rise to small corrections. In particular, for 
small spatially-isotropic perturbations around de Sitter the correction due to the terms 
quadratic in the curvature actually vanishes. These results for fields with a mass much 
larger than the inverse of the typical radius of curvature provide a specific example with 
which dark energy models based on the RG running of the cosmological constant due 
to massive fields in the decoupling regime can be compared. 

When solving the semiclassical Friedmann equation, we introduced a perturbative 
expansion in powers of (l v H) 2 . Its purpose was to obtain a fairly accurate description for 
phenomena involving length-scales much larger than the Planck length while discarding 
spurious solutions involving Planckian scales, where the effective field theory approach 
that we have been using breaks down. Furthermore, employing this kind of perturbative 
treatment made it possible to obtain explicit analytic solutions. However, truncating 
the perturbative expansion for the solutions (rather than doing so at the level of the 
equation of motion and then solving it exactly) can sometimes miss the right long- 
time behavior. In any case, our perturbative calculation would still signal when that is 
happening: one would find a solution that grows at late times up to a point where the 
perturbative expansion can no longer be applied. On the other hand, if the perturbative 
solution remains small at all times, it means that it was safe to use the perturbative 
treatment. 

An important point concerning local terms proportional to the logarithm of the 
scale factor should be emphasized. As seen in section 3, a local term proportional to the 
logarithm of the scale factor appears generically in the influence action [see equation 
(43)]. For an expanding spacetime, and especially for an exponentially expanding one 
like de Sitter, the contribution from this term can dominate at sufficiently late times over 
all the other local terms. In particular, it gives rise to a logarithmically growing term 
in the expression for the expectation value of the stress tensor which would prevent 
the existence of a self-consistent de Sitter solution and would imply an increasing 
deviation from the classical de Sitter background. However, it turns out that this 

+ It should be emphasized that a complete analysis of the backreaction problem in de Sitter spacetime 
and its stability should take into account the effect of the quantum metric fluctuations as well. Including 
the metric fluctuations is certainly a crucial aspect, and some steps in that direction are briefly discussed 
below. However, given the complexity of a completely satisfactory treatment involving the quantized 
metric perturbations, it is important to make sure that there are no significant effects even when they 
are not taken into account, especially because such effects have actually been suggested by a number 
of studies in the literature. 
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growing contribution is cancelled out by a similar contribution form the non-local term. 
This can be explicitly seen in our result for the stress tensor in the example of a massive 
field with conformal coupling. In that case the corresponding local term in the action 
is proportional to m 4 and has the same form as a cosmological constant term times the 
additional In a factor. The logarithmic term that appears explicitly in (53) is cancelled 
out by the first term on the right-hand side of (57), which contains the contributions 
from the non-local term. As a consequence, the final expression for the stress tensor, 
given by (59), exhibits no logarithmically growing term (there can be other logarithmic 
terms, but relatively suppressed by negative powers of the scale factor). Furthermore, 
we find that there is a self-consistent de Sitter solution of the backreaction equation and 
it is stable under spatially isotropic perturbations. The situation is similar in the other 
cases. For a massless field with non-conformal coupling there is a local term with the 
logarithm of the scale factor multiplying the square of the Ricci scalar.* Finally, for a 
massive field with non-conformal coupling one has the two terms of the previous two 
cases plus a term proportional to m 2 with In a multiplying the Ricci scalar. In all those 
cases there is a cancellation of the logarithmically growing contribution from the local 
terms proportional to In a and a similar contribution from the non-local term. In fact, 
this phenomenon is fairly general: one can see that something analogous also happens 
when considering spatially anisotropic and inhomogeneous perturbations [56, 65]. In 
light of these results, it may be worth revisiting related studies where the non-local 
terms where not considered. In particular, a local approximation for the effective action 
which excludes the non-local terms (but includes the local terms proportional to In a) is 
commonly employed in the literature and often leads to significant deviations from the 
de Sitter solution due to the backreaction of massive fields treated perturbatively. Our 
explicit calculation of the effective action (including the non-local terms) can contribute 
to a better understanding of the validity of the local approximation in those cases. 

In this paper we have calculated the effective action and solved the backreaction 
equations for free scalar fields. However, it is fairly straightforward to extend our 
calculations to fermions or vector fields. We expect the results to be similar to the scalar 
case, and the main conclusions to remain the same. Furthermore, it should be noted 
that our results can also be useful to compare with similar analysis of the backreaction 
due to graviton one-loop effects. Although we did not quantize the metric (and did 
not have to deal with the ambiguities associated with the gauge-fixing of the quantum 
metric perturbations), our results for the effective action have mathematically the same 
form as that in some studies of the graviton case. In particular, our influence action 
for the massless non-conformally coupled case has the same form as the effective action 

* Note that this term gives no logarithmically growing contribution to the stress tensor evaluated on 
a de Sitter background for the same reason why the local R 2 term gives a vanishing contribution to 
the stress tensor when evaluated on a de Sitter background. Nevertheless, it will give in general a 
logarithmically growing contribution when considering a background different from de Sitter (more 
specifically, whenever a background with OR ^ is considered), or when considering higher orders in 
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in [38]. Similarly, our influence action for the weakly massive and conformally coupled 
case has the same form as the dominant contribution considered in [39] if one replaces 
the cosmological constant in their result with our m 2 (their full expression including 
sub dominant terms in the infrared limit would be analogous to our weakly massive and 
non-conformally coupled case). Note that our result was derived for m 2 <C H 2 whereas 
their cosmological constant A is of the same order as H 2 . Nevertheless, this condition 
does not play any role either in the derivation of the stress tensor from the effective 
action or when solving the backreaction equation. 

We have studied the backreaction of the quantum fields on the dynamics of 
the spacetime geometry within the framework of semiclassical gravity, which can be 
understood as a mean field approximation where the mean gravitational field is described 
by a classical metric whereas its quantum fluctuations are not considered. In order to 
study the quantum fluctuations of the gravitational field one can consider the metric 
perturbations around a background geometry corresponding to the semiclassical gravity 
solution and quantize them within a low-energy effective field theory approach to 
quantum gravity [66, 67, 54, 33]. So far this approach has been mostly applied to 
weak-field problems [68], but it seems particularly interesting to extend its application 
to strong-field situations involving black holes [69, 70, 55] and cosmological spacetimes 
[71, 72]. The stochastic gravity formalism [73, 74, 52, 75, 53] can be a useful tool 
in this respect since one can prove its equivalence to a quantum treatment of the 
metric perturbations if graviton loops are neglected, which can be formally justified 
in a large N expansion for a large number of matter fields [76]. A central object in this 
formalism is the symmetrized connected two-point function of the stress tensor operator 
for the quantum matter fields, which determines the metric fluctuations induced by the 
quantum fluctuations of the matter fields. Such an object has been computed for a 
massless minimally coupled field evolving in a de Sitter background spacetime and the 
fluctuations of the stress tensor were found to be comparable to its expectation value 
[77]. Therefore, studying in detail the quantum fluctuations of the metric in this context 
constitutes a natural extension of our work worth pursuing. The results obtained here 
would still be relevant in that case because they provide the right background around 
which the metric should be perturbed and quantized. 

We close this section with a brief discussion of the relationship between our results 
and the linearization instability for metric perturbations around de Sitter spacetime 
coupled to a scalar field found in [78], where it was concluded that it is only consistent 
to consider de Sitter invariant states for the quantum field. This conclusion does not 
directly affect our analysis because we did not consider fluctuations of the metric and 
studied only the dynamics of the mean geometry, which couples to the expectation value 
of the stress tensor operator of the matter field. The expectation value of the stress 
tensor for the class of states that we have considered in this paper, which are spatially 
homogenous and isotropic, automatically satisfies the linearization stability constraint 
given by equation (44) in [78]. It is when considering the quantum fluctuations of the 
metric that the linearization stability condition imposes additional restrictions on the 
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state of the matter field because in that case the condition must be imposed on the 
n-point correlation functions of the stress tensor as well, and this implies that the state 
of the field must be de Sitter invariant. 
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Appendix 



As has been pointed out in section 3, the renormalized expectation value of the trace of 
the stress tensor may suffer from initial time divergences unless we choose an appropriate 
state for the quantum field. In this appendix we will carefully examine this point. For 
simplicity, we will restrict ourselves to the massless case, m = 0, v ^ 0. However, no 
assumptions on the size of the v parameter will be made here, so the following results 
apply for both the weakly and the strongly non-conformal regimes. Throughout the 
appendix we will follow closely the definitions and results for fourth-order adiabatic 
vacua of [46]. 

We will make use of an alternative expression for the renormalized expectation 
value of the trace of the stress tensor. Consider the mode functions solution of (7) 
with any initial conditions consistent with the Wronskian condition (8), not necessarily 
of the form (10). The bare expectation value of the trace of the stress tensor in the 
vacuum state associated with this set of modes can be written as (see, for instance, [79]) 



(TP 



6is 



2n 2 a 4 



dkk 2 



h\f k \ 2 + 2h\f k 



d\f k 
drj 



\fk\ 2 + ^l\fk\ 



(Al) 



where we have used the notations Q 2 , = k 2 + M 2 and h = a/a. In order to renormalize 
this expectation value, it is convenient to introduce the WKB mode functions. The 
most general mode function satisfying the Wronskian condition (8) can be written as 



fk(v) = 



Vim 



exp 



drfW h {rf) 



Vi 



(A.2) 



with real and positive time-dependent frequency W k . [We are assuming that f\ is real 
and positive at the initial time, which implies no loss of generality because of global 
phase freedom.] With this change of variables, the mode equation (7) takes the form 

Wk_3(W h 1 
W k 2\W kf 

Solving this equation iteratively, starting with W k = ^4, one obtains an adiabatic 
expansion for W k . This is the WKB solution of (A. 3), and the mode associated with it 



W 2 k 



k 2 



(A.3) 
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through (A. 2) is a WKB mode. The WKB frequency admits the following expansion in 
inverse powers of k: 



W k = k + 



M 2 M 4 + (d 2 M 2 /dri 2 ) 



+ o(k- 



(A.4) 



2k 8k 3 

The adiabatic approximation of fourth order, which contains up to fourth-order 
derivatives of the scale factor, is obtained by truncating the iterative procedure after 
the second iteration, 

1 £l k 3 / 

2 Qk 2 \ Qk 



W, 



(4) 



(A.5) 



and differs from the exact solution Wk(rf) only in 0(k~ 5 ) terms. The bare expectation 
value (A.l) can now be renormalized by the so-called adiabatic subtraction procedure 
[80, 81, 82], which consists of subtracting the same expectation value computed with 
the fourth-order adiabatic approximation to the WKB modes: 

{TPren={T»)-{T^\ (A.6) 

where (T^ 1 )^ is the trace (A.l) associated with the modes (A. 2), with time-dependent 
frequency W k given by (A.5). This procedure has been shown to be equivalent to 
covariant renormalization methods specialized to RW spacetimes [83, 84]. 



Physical vacua 

What kind of initial conditions should the modes {f k } satisfy in order for the 
renormalized expectation value of the trace (T^) ren to be finite at the initial time 
?7i? According to (A. 2), the most general initial conditions that fulfill the Wronskian 
condition (8) are 

fk(m) 



\/2w k 

fk(m) = (~iw k + y) fk(Vi), ( A - 7 ) 

where w k and v k are two arbitrary real functions of k (with w k > 0). Substituting these 
equalities into (A.l) and (A.6), we get that the renormalized expectation value of the 
trace at the initial time is 
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where the subscript "i" indicates that the function is evaluated at the initial time rji. 
This integral is finite provided that the integrand falls off faster than k~ x as k — > oo, 
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which is achieved if 

\w k -W ki \ <0(AT 3 ) 



v k + 



< 0(k~ 2 ) 



(A.9) 



when k — > oo. Here we have used the fact that W k and wjf^ differ only by 0(k~ 5 ) terms. 
These are the constraints that select, among all the possible homogeneous and isotropic 
vacua, the physical ones, namely the ones that keep the trace finite at the initial time. 
They correspond to the so-called fourth-order adiabatic vacua. 



Our vacua 

Let us now see how these constraints affect our particular class of vacuum states. In 
other words: among our vacua, which are the physical ones? Our vacuum states are 
associated with the initial conditions defined by (9)-(12). Note that, in fact, in the limit 
k — > oo the modes {f k } are just the WKB modes of equation (11), 



1 



exp 



dr/W^rf) 
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where is given by (A. 4) with Mq, replacing M. The reason for this is that, for 
k — > oo, both sides of the equation above satisfy the mode equation (11), and both 
sides reduce to standard plane waves when 77 — > — 00. Therefore, our particular initial 
conditions can be written in the form (A. 7), with w k and v k given in this case by 



v k = 



(A.H) 
(A.12) 



00. The constraints (A.9), together with the expansion (A. 4), then imply 

< 0(AT 3 ) 
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^[M 2 , - M 2 } - JL [ M |, - Mt + {d 2 Ml/d V 2 ) i - {d 2 M 2 /drf).] + 0(AT 5 ) 
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2k 2 
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< 0(k~ 2 ). 



(A.13) 



In order for these inequalities to be satisfied, we must impose 

Ml, = M 2 (dMl/drj). = (dM 2 /d V ). {d 2 Ml/dr 1 2 ). = (d 2 M 2 /d V 2 ). . 

which imply precisely the conditions in (70), which were obtained by other means in 
section 3. Therefore, the requirement that the scale factors and a are matched 
smoothly enough selects, among all the states considered in this paper, the fourth-order 
adiabatic vacua. 

We close this appendix by mentioning an alternative method of calculating the 
influence action (but entirely equivalent to that of section 3) provided in [85]. The 
approach, which is based on decomposing the field in spatial Fourier modes, computing 
the unitary evolution operator for each mode perturbatively in the interaction picture, 
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and summing over all the modes at the end, can be useful when considering more general 
initial states at a finite initial time ^ which are not necessarily of the form described in 
section 2.1. 
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